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ABSTRACT 



We present extensive tables of the Atomic Partition Function for iron ions, Fe I - 
Fe X, and discuss details of the computational method. Partition functions are given 
in wide range of temperatures, 10 3 K < T < 10 6 K, and lowering of ionization energy 
(0.001 eV < LIE < 5.0 eV). Our APF take into account all energy levels predicted by 
quantum mechanics, including autoionization levels. The tables can be applied for the 
computations of model stellar atmospheres and theoretical spectra over all existing spectral 
and luminosity classes. 

Our tables are available at |http : / / www . astrouw . edu . pl/^acta/ acta . html] (Acta 



Astronomica Archive), and [http : //draco .uni . opole .pl/Halenka. html . 
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1. Introduction 

Iron is the most abundant and important heavy element in astrophysical plasmas. Due 
to its relatively high abundance in solar type composition and low ionization energy, iron 
ions are important contributors of free electrons in stellar atmospheres of population II 
at low and moderate effective temperatures (spectral types F and lower). In atmospheres 
of these types hydrogen and helium atoms remain almost perfectly neutral, and therefore 
iron donors of free electrons can be critically important for the populations and opacity of 
negative ions H~ and He", which in turn determine the whole atmospheric structure and 
visual/infrared spectra of these most numerous classes of stars. 

Also on the opposite end of the Hertzsprung-Russell diagram, various ions of iron 
contribute to the stratification of monochromatic bound-free and line opacities in stellar 
atmospheres of the highest effective temperatures, including both main sequence and white 
dwarf stars. Highly ionized iron ions efficiently absorb extreme ultraviolet and X-ray 
photons, where the outgoing flux of radiation is most intense in these objects. Therefore 
iron plays a crucial role in determining the structure of atmospheres and appearance of 
the radiation spectra of the hottest stars. 

Both the above examples show the importance of construction of the partition functions 
for various iron ions, which are internally homogeneous and can determine exactly 
populations of iron ions over very wide range of temperatures and electron concentrations. 
Our paper fulfils such a demand. 

Values of the partition function U for a particular ion depend on gas temperature T, 
and the local electron concentration N e . Temperature enters partition function by its 
definition (Griem 1964, Drawin and Felenbok 1965, Traving et al. 1966) 

U(T,N e ) = ^g^M-Ei/kT) , (1) 

i=l 

where the sum is taken over discrete energy levels of statistical weights i and excitation 
energies E { . In vacuum (N e = 0) the number of bound levels i max is infinite, and therefore 
the above series diverges. 

In real plasma, however, interaction between the atom of interest and surrounding 
free electrons and ions cause, that bound levels of very high excitation energies move to 
continuum and no longer contribute to the partition function. Therefore the series in Eq. 
(1) reduces to finite number of terms, and value of U is also finite and is strongly dependent 
on the electron concentration N e . In general, the larger is N e the lower is both number of 
bound energy levels i max in Eq. (1), and the value of U. 

There exist a large number of papers, which present tables of rather approximate (or 
even schematic) partition functions (APF) for elements, including few iron ions (Drawin 
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and Felenbok 1965; Traving et al. 1966; Irwin 1981, for example). The latter paper presents 
fitting formulae for APF of Fe I - Fe III, for temperatures T < 16000 K. The widely used 
computer code Tlusty 195 for computations of NLTE model stellar atmospheres (Hubeny 
and Lanz 1992, 1995) contains FORTRAN subroutine computing APF of Fe IV - Fe IX 
by direct summation over all energy levels observed in these ions in 1971. Therefore the 
set of energy levels was not complete in TLUSTY195. 

Quality of given tables of partition functions depends on (i) accuracy and completeness 
of energy levels included in Eq. (1), and (ii) realiability of the assumed theory of emitter- 
plasma interactions. Emitter-plasma interactions cause that the series in Eq. (1) is finite. 
Unfortunately, none of the currently existing theories describe correctly effects of charged 
particles in plasma on the atomic partition functions, cf. also Hummer and Mihalas (1988). 

We have computed and present tables of the partition functions taking into account 
all energy levels predicted by quantum mechanics, including also levels lying above the so 
called normal ionization energy (autoionization levels). Moreover, our computations are 
more physically correct, since values of APF depend on both temperature T and electron 
concentration N e . The latter is related to the lowering of ionization energy (LIE). The 
method of APF calculations is described in the following Section. This method was also 
presented in series of papers by Halenka and Grabowski (1977, 1984, 1986), Halenka (1988, 
1989), and Madej et al. (1999). Paper by Halenka et al. (2001, hereafter Paper I) presented 
already extensive tables of APF for nickel, Ni I - Ni X, which were computed in the same 
way. 



2. Computations of the atomic partition functions 

We have appended contribution to the Atomic Partition Function given by autoionizing 
levels. We replace the classical definition of APF, Eq. (1), by more realistic expression 
introduced by Halenka and Grabowski (1977) as follows 

Pmax 

Hp) max Pmax 

U^(T,N e ) = E E g^expi-E^/kT) = £ U^(T,N e ) , (2) 

p—1 i—1 p—1 

Here the set (pi) with the indices p and i (ordering of levels from the ground level upwards 
in energy scale) describes an eigenstate of the atom in the r-th ionization state. Index 
i represents three quantum numbers (nlj) of the optical electron, and p represents the 
quantum state of the atomic core. Index i(p) max is the number of all bound energy levels, 
g$ and E$ denote statistical weight and excitation energy of the i-th state, in the sequence 
based on the p-th parent level. Numbers i(p) max result from the inequality 



E%> < Eg - AE^ , 



(3) 
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where Ai? (r) denotes the lowering of the ionization energy LIE, and Ej^ is the ionization 
energy in the p-th level sequence. The latter quantity is equal to the sum 

E (r) = E (r) + £(r+l) 

poo loo 1 p V / 

The quantity E^ +1) denotes the energy of the atomic core after ionization, r — > r + 1. 
Index p TOaa . is the number of different parent levels which can be realized in given physical 
conditions. The number p max results from the inequality similar to Eq. (3), written for 
the (r + l)-th ionization state. Since for a fixed value of p the number k is assigned 
unambiguously, then the upper limit of E^ for the fc-fold excitation can be written as 
follows. Numbers k = 1,2, . . . , Z — r, where Z denotes the atomic number 

E^<tE[T- 1} - (5) 

s=l 

The above definition of APF requires, that for a fixed ionization state r the initial set of 
observed and calculated energy levels must be split into level sequences according to their 
attachment to different parent terms. Our initial set of levels was taken from catalogue of 
energy levels of iron ions from the CD-ROM by Kurucz (1994), which contained both the 
'observed' and computed energy levels (of the order 1000 for each ionization state r). We 
have also added to our APF's contribution from many energy levels predicted by quantum 
mechanics, which were missed in CD-ROM by Kurucz. Exact description of the method 
of adding of missing levels was given by Halenka and Grabowski (1977, 1984). 

Given term of a configuration with equivalent electrons usually belongs simultaneously 
to few level sequences defined by different parent terms. To calculate the contribution of a 
fixed term i to the partial APF, U^ r) in Eq. (2), one must multiply the statistical weights 
of the z-th term by their coefficients of the fractional parentage G qi , which were normalized 
to 1. We have assigned coefficients of the fractional parentage according to the well known 
procedure (cf. Slater 1960, vol. II; Sobelman 1979): 

G pi = G a > L i s i. aLS = (aL's'; aLSf/ ^ ' a L' S' ; aLSf. (6) 

For a fixed ionization degree r, all level sequences were divided into three groups: a {r \ 
/3 (r) , and 7 (r) . 

Group consists of level sequences numbered by p = 1,2, ...pi, where all energy 
levels (energies and statistical weights) are known. Usually the total number of such level 
sequences does not exceed 20. In most cases these level sequences belong to the ground 
parent configuration. Their contribution to the APF is given by summation of terms in 
Eq. (2). 

Group consists of level sequences, in which only part of energy levels are known, 
from i = 1 up i\(p)- Number of such level sequences, p = p 1 + 1, . . . ,p 2 , never exceeds one 
hundred. Contributions of /9 (r) level sequences to the APF is given by summation up to 
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level ii(p). Remaining levels are included approximately, assuming their similarity to some 
of level sequence of the group a (r) , so called reference level sequence. Thus, for pi < p < p 2 , 
values of U { p r) can be calculated according to the relation: 

hip) 

E u^= £ E^ ) «p(- s £V*t)+ E ACtf°, (7) 

pEterm s pGterm s i—1 pGterm s 



where 



V An (r) ~ V Ar r(r) peterm s rri 

^ p ~ ^ p V r7(r+i)pW-R (r+1) /Fn ' 1 J 

term s pGterm q 2^ i/p ; «Xp^ £/ pi / ft-J. J 



pet. 

pGterm q 



and 



Here (/ (<? < Pi) represents a reference level sequence. 

Group 7 (r) is the set of level sequences, in which no levels are known. Their share in 
the APF can be estimated from the expression 



Pmax Z—/ p 

x 



\ ^ Jj{r) peterm q 



£1 P E 9^eM-E£ +1) /kT) 

pGterm q 



E 

^ (r+1) -E^ r+1)ex p(-^ +1) /^) 

p=l 



(8) 



The relation Eq. (7') is strictly valid if, at a fixed quantum number J of the optical electron, 
change s — > (/ of the state of the atomic core does not alter the electron-core interaction 
energy. In such a case energy levels of the s sequence are shifted as a whole relatively 
to q sequence by energy |i?( r+1 ) — E^ +v> \, and the statistical weighs satisfy the relation 

9^ / g q r) = g [ ; +1) / g [ ; +1) ■ 

The above condition is never strictly realized. However, the error introduced into 
by the Eq. (7') reaches minimum when the maximum bound energies in both considered 
level sequences (s and q) reach similar values. It should be noted that values of £7 (r) are 
of physical interest only in such range of temperatures T, when E^/kT 3> 1. Since for 
p > pi we have E$ > E^ (except for few of the lowest levels) , then the error contributed 
by AU$ r) to U (e) in Eq. (7') is always small. 



3. Results 



We have computed APF for Fe I - Fe X ions according to above procedure. Exemplary 
results shown in Fig. 1, which displays run of our APF for Fe IV and temperatures 
corresponding to atmospheres of hot white dwarf stars, for various values of LIE taken 
as free parameter (solid lines). Numerical values of log APF for Fe IV are also listed in 
Table 1. In Fig. 1 we have appended run of partition functions computed on the basis of 
the observed levels only (single dashed line). One can easily note, that our recommended 
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partition functions are significantly larger than the latter functions, and apparently tend 
to diverge at any fixed T and iV e — > 0. 

Partition functions for iron are arranged in 10 ASCII tables, where each table 
corresponds to a single ion, which is identical with the organisation of Ni I - Ni X partition 
functions of Paper I. Full set of these data is available from Acta Astronomica Archive, 



http : / / www . astrouw . edu . pi/ ^acta/ acta . html 



or [http : / / draco . uni . opole . pl/Halenka . html . 

Entries of tables are decimal logarithms of the APF. They are tabulated at 62 discrete 
temperatures, spaced at nonequidistant intervals, 10 3 <T< 10 6 K, and at 9 arbitrarily 
assumed values of the lowering of ionization energy {LIE = 0.001, 0.003, 0.010, 0.030, 
0.100, 0.300, 1.000, 3.000, and 5.000 eV). Both T and LIE points remain identical in all 
10 tables of APF, to ensure homogeneity of the data and compatibility with Paper I. 

Very wide range of both parameters T and LIE assumed in our tables ensures, that they 
cover practically all the conditions expected in stellar atmospheres of any type, reaching 
from the coldest dwarfs up to the hottest accreting white dwarf stars. 

Tables discussed in this paper are based on the most complete set of energy levels 
actually available for iron ions. Moreover, our partition functions are sensitive to plasma 
interactions, since they strongly depend on the lowering of ionization energy. Our APF 
tend to diverge for the lowering of ionization energy approaching zero, and are perfectly 
compatible with the APF for nickel ions (cf. Paper I). 
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log(T) (K) 

Fig. 1 : Run of partition functions of Fe IV as function of gas 
temperature T and various parameters LIE. Solid lines represent our 
recommended results, whereas dashed line represents partition function 
computed from the observed levels only. 
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Table 1 

Atomic partition functions for Fc IV (decimal logarithms) 
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